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Abstract 

We investigate structure functions in the 2-dimensional (asymptotically free) non- 
linear O(ra) cr— models using the non-perturbative S-matrix bootstrap program. In 
particular the exact small (Bjorken) x behavior is derived. Structure functions in the 
special case of the n = 3 model are accurately computed over the whole x range for 
-q 2 /M 2 < 10 5 , and some moments are compared with results from renormalized 
perturbation theory. Some results concerning the structure functions in the 1/n 
approximation are also presented. 



1 Introduction 



In this paper we study structure functions in the asymptotically free O(n) sigma 
models in two dimensions. Due to the integrability of the model one has powerful 
tools to study various non-perturbative properties. In particular one can derive the 
exact small x behavior (for all q 2 ) and for the case of n = 3 compute structure 
functions precisely up to very large values of q 2 . Despite the fact that there are 
no transverse directions, the structure functions have a rather rich and non-trivial 
behavior. In a previous letter [Q we summarized our results and speculated on the 
possibility of discovering some similar structural features in QCD. 

The purpose of this paper is to supply the derivation of the results presented 
in pQ. This paper is organized as follows. In the next section we give some basic 
definitions of the correlation functions of interest. In sect. 3 we give the derivation of 
the (rather universal) exact small x behavior. Section 4 deals with certain general 
aspects concerning the relation of the high q 2 behavior of moments of structure 
functions to the operator product expansion (OPE). In sect. 5 we consider the OPE 
for the cases of two spin fields and two currents in the framework of perturbation 
theory. More detailed results on the structure functions for the special of n = 3 are 
presented in sect. 6. Finally in sect. 7 we consider computations in the leading order 
of the 1/n expansion. Many technicalities and some conventions can be found in the 
appendices. 

2 0(n) model and structure functions 

The O(n) cr-model in 2d (formally described by the Lagrangian (|5.1j0 is pertur- 
batively asymptotically free for n > 3. A special property is that these models 
have an infinite number of local |2j and non-local [Hj classical conservation laws 
which survive quantization. At the quantum level they imply absence of particle 
production. Assuming the spectrum to consist of one stable 0(n)-vector multiplet 
of mass M, the S-matrix has been proposed long ago by the Zamolodchikovs [3]. 
Form factors of local operators can be computed using general principles j r )l(ij . The 
S-matrix bootstrap program for the construction of correlation functions involves 
summing the contributions over all intermediate states [7j. The possible equivalence 
of this construction to the continuum limit of the lattice regularized theory has been 
investigated in ref. |Hj. 
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2.1 Current and spin operators, 2 point functions 

The normalization of the conserved 0(n) current operator J^ b (x) (a, b = 1, ... n) is 
fixed e.g. by the equal time commutation relation with the spin field <l? c (y): 

Jg b (0,x l ),$ c (0,y 1 )] = it? d 5(x 1 -y l )$ d (0,y l ), (2.1) 

where the matrices t ab given in 1)A,2|) yield the vector representation of the O(n) Lie 
algebra. Its matrix elements are 1 

{0\jf (0)1^,0!;. ..;a r ,e r ) = -ie^P? f^fa,. . . ,9 r ) . (2.2) 

Here the number of particles, r, has to be even and the form factors depend 
on the rapidity differences only, making Lorentz invariance and current conservation 
manifest. The normalization of the r-particle states, the corresponding completeness 
relations, and other undefined kinematics encountered are given in Appendix B. We 
define 

E f£Mh,---,0r)f5Ufa,...,9r) = (^S*-S°f8*>) (2.3) 

ai ...a r 

where C^ r \u) is a symmetric function of the rapidity differences. 

The normalization of the spin operator $ a (x) is fixed by its 1— particle matrix 
element: 

(0|$ a (0)|6,#) = 8 ab . (2.4) 

Its r-particle matrix elements (r odd) are defined by 

<O|* B (O)|oi,0i; . . .;a r ,9 r ) = K f a ai ... a M, ...,8 r ), (2.5) 

where the form factors , depend on the rapidity differences only and the overall 
factor A n is defined for later convenience. We choose 

A 3 = 4=> A n = l (n>3). (2.6) 

V7T 



The analog of (|2.3|) for odd r is 

E • • • A) fataM, ...,6 r )= 5 ab C«(«) . (2.7) 
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We now make some further definitions: 



JW W= CW W „ (2, 



and for s = 0, 1, 



I s (z) = J2l^ 2k+1+s \z), (2.10) 



fc=0 



00 ^ 
A s (*) = 53A( 2fc+1+ -)(z) = -z 2 gl s (z). (2.11) 

fe=0 

The invariant functions I s are related to the 2-point functions of the current 
and spin field operators by |Hj 

(0\T*J?(xW(y)\0) = 

valid up to contact terms and 

(0\T<S> a (x)$ b (y)\0) = A 2 n 8 ab J ^^ x -y\-i)h{-p 2 - is) . (2.13) 



For n = 3 we can define J„{x) = ^e afec J^ c (x) and instead of 1)2.2)1 we have 



In 1)2.12(1 T* denotes the covariantized T-product. 

T«(x) = \e abc J b < 

(0|J2(0)|oi,ei; . . . ;a r ,9 r ) = -ie^ a P r a /^...^i, . . . A) • (2.14) 
In this case instead of 1)2. 3|) we can use p. 7)1 also for r even. 

2.2 Structure functions, moments 

The central object in DIS theory is 

W-^ ;cde/ (p,g) =vr ^(a,p|J^(0)|r)(r|J^(0)|6,p)^ 2 )(p + g-P r ), (2.15) 
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where q 2 < 0. We will use the parameterization 

q 2 = -4k 2 M 2 (2.16) 

and the Bjorken variable 

x = f— . (2.17) 

2(pq) 1 ; 

Using Lorentz and O(n) invariance we have 

W^ f ip, q) = Uu " ^) E Rf' Cd£f Mq 2 ,x) , (2.18) 
^ ^ ' 1=0 

where the projectors R[ corresponding to the 3 invariant t-channel "isospins" are 
defined in Appendix A. Note that in 2 dimensions there is only one independent 
structure function for each isospin channel. 

Similarly we define the structure functions corresponding to the spin operator 
through 

^ cd (p,q) =-V Y,(^ C mr)(r\^ d (0)\b,p)5^(p + q-P r ), (2.19) 

r 

and 

2 

E a6 ^(p,g) = A 2 n £ P^ cd m(q 2 ,x) , (2.20) 

1=0 

where the t-channel projectors Pi for the vector representation are given in (|A.5|1 — 

(EH) . 

Separating the r-particle contributions we have 

wi(q 2 ,x)= E w t\l 2 ^ x ) and wi(q 2 ,x)= E ^ (r) (g 2 ,x) (2.21) 

r odd r even 

with 

w^(q 2 ,x) = ^ F |°°dA | T>u^5^{p + q-P r )j\ r \e). (2.22) 
Here p = (M, 0) and for r odd 

J / r) W = ^ E < ;Cde/ /at...a r (^,0i,...,e r )/£ ar (ivr^ 1 ^ (2.23) 

abcdef 

a\ ...a r 
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while for r even 

J i T \°) = ^ E Pt b ' Cd fL 1 ^^^---^r)ftt..a r (^0 u ... : 9 r ). (2.24) 



abed 

ai ...a r 



By doing the A-integration we can further simplify (|2.22[) : 



wf\q 2 ,x) 



2k 2 



(4tt)' 



r-l 



where 



b 



4-i 



4k 2 



J, (r) (A,...,A0, (2.25) 



X f X ^ 

^ + 4^/- ln l 1 - X+ 4^) 



We define the structure function moments by 
-1 

JV-l 



M Z;iV (^) = / dxx^" 1 ^ 2 ,*) and M Z;7V (c/ 2 ) = / dx x 7 ^ 1 ^(g 2 , x) 
and similarly for fixed particle number 

M^(q 2 ) = C dxx N - 1 w i [\q\x). 



(2.26) 
(2.27) 

(2.28) 
(2.29) 



Obviously 

M, N (Q 2 ) = E M S(? 2 ) 



and M^(g 2 ) = ^(g 2 ) . (2.30) 



r odd 



The r-particle moments can also be calculated directly from (|2.25|) : 



O 2 ) 



1 



2(4vr)'- 1 
4k 2 



X N+1 4 r \P!,...^r) 



4k 2 + $ - 1 ' 
For n = 3 Q2.15[) can be written as 



(2.31) 
(2.32) 

(2.33) 



and (|2.18|) becomes 



w$*m = ( v,u - ^ ) E JfWv 



(2.34) 



In this case (j2.24|) is valid for odd as well as even r values (with tti = 21 + 1) 
The 2-particle form factor can be written 



i — c/2 J \6 ac S M -S ad 6 bc 



with 



l anh - ex]") { ~ - I — 



l-e 



l + e* 



sin 2 ([i7r - 6}t/2ir) 
sinhi 



The 1-particle contribution to the structure functions is then given by 



w\ ; (q 2 ,x) = mi 8{x — 1) \4>{m — a)\' 



where 



sinh — = k . 
2 



and 



mo = 1, mi 



-m 2 



(2.35) 

(2.36) 

(2.37) 
(2.38) 



3 2c? structure functions at small x 



In this section we derive a general formula describing the asymptotic behavior of the 
O(ra) model structure functions at small x values. The derivation is based on general 
properties of the form factors and the scattering matrix elements and therefore the 
behavior we find here is expected to hold in other 2d integrable models as well. 



For small x — > the variable b in (|2,27|) behaves as 

b = x + 0(x 2 ), (3.1) 
and if we do the u r -i-fntegration in (j2.25|) with the help of the delta function we get 

u r _! = In - 2v { l~ l) + O(x) , (3.2) 

and further 

TTI/. 2 , _ r)// 2 zIk- 2 

(3.3) 



(r) %H r -i „ , 2\ 



Ak 2 

Now putting all the above together we have 

Sr)(„ 2 - - • x 



da 2 Ak 2 

^^ = ±L + o(i) 

OUr-l X 



w} >(q ,x) 



where 



2(4vr) 



r-l 



JY'i-e, -p-e + Pr-i, ...,-P-e + Pi), (3.4) 



l + f^)x + Q(x 2 ) 



and 



/? = lnx + 0(l) (3.5) 
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and Pi are the variables for r — 1 particles. 

Eqns. and are both of the form 

jf\0)= Y, ^^/« 1 ..-,(^«l»---^r)/&..^,^...,flr), (3-6) 
ai...ay 

where for the current case 

c ab;AB = J_ R ab;cdef ^ A ~ cd , B ~ e/ , (3.7) 



and for the spin case 



C flMB = — R a6;cd with A ~ c , B ~ d . (3.8) 

717 



With this notation we can write Jj in (|3.4() as 



abAB 



(3.9) 



The crucial point now is that since e is small and /3 is large we can here use (|L).4|) 
which follows from general principles encoded in the Smirnov axioms. In leading 
order we get 

T (r) 

(47T) 2 

\ A /~iab\AB j.aa r j.ba r fA' /a a \ 

J l — ( n -2) 2 e 2 8 2 AA' BB' J a r -\...ai \Pl 1' ■ • • i PI) 

V ' P abAB (3.10) 

a\...a r v / 

which can be further simplified with the help of (f%~20l) . (fQ2jl . (tO|l and tf2~7^ 

leading to 

where the constants G/ are equal to V/ and T] for the spin and current cases respec- 
tively given in (TOP) and (TQ3l) . and further 



u> ' V^) =- t— ' ^ / , (3.12) 

xln 2 x (4^-3(7! - 2)2 J (4« 2 + M?-i) 

which can also be written as 

«fV,*) ^T^^H 2 ), (3.13) 
x In x (n - 2)^ 
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where the Adler functions A^' were defined in (]2.9|) , 

The final results for the complete structure functions are 



w t {q 2 ,x) 



1 2ttT, 




(3.14) 



xln 2 x (n - 2) 2 



and 




1 27rV^ 




(3.15) 



xln 2 x {n - 2) 2 



Note that the structure of the asymptotic small x behavior, factorizing a part charac- 
teristic to the target and a part described by the vacuum 2-point function, is rather 
universal being independent of the operator, independent of n, and independent of 
the isospin channel. 

This completes the derivation of the exact small x asymptotics first announced 
in ref. pQ. The question of possible lessons that can be learned for QCD was ad- 
dressed in the latter reference and will not be repeated here. 

4 The operator product expansion 

In the O(n) cr-models there does not seem to be a simple parton picture. This is 
even so for the case n = 3 where the model is equivalent to the CP 1 model. For 
although this model is formulated in terms of a complex doublet of fields which 
are analogous to quarks in that they are confined, it seems that they do not play 
a role more similar to partons than the elementary bare spin fields in the original 
formulation 2 . The question is related to that of understanding what are (if any) 
the "ultra-particles" in the sense of Buchholz and Verch [10] . or to the associated 
question as to whether the <7-models have an underlying conformal field theory. 
Although an intuitive parton description with suggestive DGLAP equations 



(where pi(z,q 2 ) would be the corresponding splitting functions) is still missing in 
these models, we still have the machinery of the operator product expansion (OPE) 
to give us information on the evolution of the moments ()2.28j) at large — q 2 . 

The OPE in the sigma model is surprisingly involved and hence we have decided 
to present the material as follows. In the next subsection we first summarize the 
results; readers who would prefer to skip the derivations can then jump to sect. 6. 

2 Perhaps the peculiar threshold behavior discussed in sect. 6 is explained by the fact that (as 
opposed to QCD) with some probability the O(n) particle can consist of a single point-like parton 
that carries the same quantum numbers. 




wi(q 2 ,x)= — Pl{x/y,q 2 )wi(q 2 ,y) 




(4.1) 
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The general structure of the product of two local operators (in this case the spins and 
currents) is described in the remaining part of this section. Our analysis extends that 
initiated e.g. in refs. So far too little is rigorously known about the detailed 

structure of the OPE from the general principles of the bootstrap approach to obtain 
the explicit results below. The extra required information is however supplied in the 
framework of renormalized perturbation theory which is presented in sect. 5. Some 
comparisons of the moments with those from the bootstrap approach at high —q 2 
are presented in subsect. 6.4. 

4.1 Summary of results on the moments 

For the current (N even) moments in the isospin channel we have 

M , N {q 2 ) = W,, Nl ^ T) |l + -l^A(g 2 )+0(A 2 )| , iV>2, (4.2) 
where X(q 2 ) is an effective running coupling function defined through 

' ' ' lnA(g 2 )=ln^3, (4.3) 



X(q 2 ) n-2 A- MS 

and the Wq-n are renormalization group invariant, non-perturbative constants, cor- 
responding to the matrix elements of spin N operators. In the N = 2 case this is the 
energy-momentum tensor operator T^ u for which we know the constant explicitly 

(a,p\T^(0)\b,p) = W . 2 p^5 ab , W Q , 2 = 2. (4.4) 

In particular the "momentum sum rule" follows: 

n — 2 

M 0;2 (-oo) = . (4.5) 

n — \ 

Note that all the isospin moments tend to constants as — q 2 — > oo. As a conse- 
quence these current structure functions in the O(n) models obey Bjorken scaling. 
Computations in the n = 3 model, (see sect. 6 and in particular Fig. 3), indicate that 
the resulting limiting scaling functions are non-trivial. This is a special property of 
these models and we conjecture that this is due to the existence of an infinite set of 
local conserved quantities 

In the isospin I = 1 channel for odd moments N > 3 we can only say that 

Mi-, N (q 2 ) = Wi-n \{q 2 )^ + . . . , N > 3 , (4.6) 
but in the special case TV = 1 we have 

M l , l (q 2 )= l - |l-_l_A(g 2 )+0(A 2 )} , (4.7) 
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where the constant is known through the current normalization 



(a,p\jf(0)\b,p) = -4ip,Pf> cd . (4.8) 
From this follows the analogy to the Adler sum rule in QCD: 

Afi ; i(-oo) = i . (4.9) 
For the spin field isospin moments we have 

MoMl 2 ) = ^ 7r2 "f" A ( g 2 )^{l + 0(A)}, n>4, (4.10) 
(n — 2) z 

MoMQ 2 ) = ^{1 + A( (7 2 ) + 0(A 2 )} , n = 3, (4.11) 

where the non-perturbative constants Wq-n are the same as for the current, and 
where C n is the non-perturbative constant appearing in the short distance expansion 

(0|<l> a (y)<t> b (0)|0> ~ Cn 5 ab (-lnM\y\)^ . (4.12) 

So far the value of C n is not known for general n; for the case n = 3 a (well tested) 
conjecture based on scaling j^j gives 

Ck = (4.13) 

and we know for n = oo: 

c ~ = h- (4 - 14) 

We see that only for the case n = 3 do the moments of the field / = structure 
function have the same leading asymptotic behavior as those of the current. 
For the isospin I = 1 field (odd) moments we find to leading order PT 

M 1;1 (q 2 ) = M 0;2 ((Z 2 ), (4.15) 

M 1;JV ( g 2 ) = W 1]N X(q 2 )^ {l + O (A^)} , N > 3, (4.16) 

where there is in general no obvious relation between the Wi-n and the constants 
occurring in (|4.6|) . except for n = 3 where they are equal (W\-n = Wi-n > n = 3). 
For isospin / = 2 moments we obtain (for all n > 3): 

M 2 . N (q 2 ) = W 2 -nK<1 2 ) 2 {l + O (A^) } . (4.17) 
Finally using the exact ratio of the mass to the A-parameter 

JW_ (8/e)V(- 2 ) 

A^s r[l + l/(n - 2)] ' 1 ' 1 

obtained by Hasenfratz, Maggiore and Niedermayer [T5|, the perturbative results 
can be plotted as functions of —q 2 /M 2 . 
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4.2 Dispersion relations 

For the discussion of the OPE it is convenient to work in the Euclidean formalism. 
For local operators A we have 

A(x°, x 1 ) = e ^ Hx °- p ^ A(0, 0) e-*^ -^ 1 ) , (4.19) 

where H, P are the time and spatial translation operators. We can similarly define 
Euclidean translation by 

A E {y u y 2 ) = e Hy2 - iPyi A(0, 0) e' Hy2+iPyi , (4.20) 

which is formally A(—iy2,yi)- For Euclidean vectors (and similarly for tensors) we 
define 

V 2 = -iV . (4.21) 
The Euclidean time ordering is defined as 



T E (A E (y 1 ,y 2 )B E (z 1 ,z 2 )) = 0(y 2 - z 2 ) A E (yi, y 2 )B E (z 1 , z 2 ) 

+ Q(z2-y2)B E (z 1 ,z 2 )A E (y 1 ,y 2 ), 



(4.22) 



and the connected part of the product of two operators is 

(O l 2 ) c = O x 2 - <0|OiO 2 |0) - Oi|0)(0|O 2 . (4.23) 
We now define Euclidean functions for the currents: 



\ J dy 1 dy 2 e^ + ^)(a,0|^(j^ E (y 1 ,y 2 )J^(0,0))j&,0) 

2 

= (QfiQu - Q\u) Rf cdef n{Q\Q2) , 



(4.24) 



where stands for covariantized Euclidean time ordering, i.e. some non-covariant 
terms proportional to delta functions of the Euclidean time difference (and deriva- 
tives of the delta function) are dropped. 
Similarly for the spin field: 

\ j d yi dy 2 e^ 1+ ^ 2 )(a,0|r E ($ cE ( yi , y 2 ) <K d (0, 0)) J6,0) 

2 (4.25) 
= A 2 J2P? b ' cd f l (Q 2 ,Q 2 ). 

1=0 
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The T\ and 77, as functions of Q2 at fixed real Q 2 are real analytic 

ti(Q 2 ,Q2)* = n{Q\-Q* 2 ), f,(Q 2 ,Q 2 )* = r l (Q 2 ,-Q* 2 ), (4.26) 

and obey the crossing properties 

ti(Q 2 ,Q 2 ) = (-l)S(Q 2 ,-Q 2 ), f l (Q 2 ,Q 2 ) = (-l) l f l (Q 2 ,-Q 2 ). (4.27) 

Further they have cuts along parts of the imaginary axis (with poles at Q2 = 
±iQ 2 /2M), and the discontinuities across the cuts are simply related to the structure 
functions: 

wi(-Q 2 ,x) = 9li mTl (Q2 :£ _i_^L\ ( 4. 28 ) 



vr V 2Mx 

wi{-Q 2 ,x) = 9-lmfi (Q 2 ,e-i^- ) . (4.29) 



vr 'V 2Mx 

Concerning the general singularity structure in the complex Q2 plane away from 
the imaginary axis, little more is rigorously known except that the current function 
7~i has poles on the real Q2 axis originating from the contribution from 1-particle 
intermediate states. The contribution from 1-particle states is easily computed: 

7 i l ^(0 9 ,ft) = - 5S ^^-*Wir + fc)x 

1 ( _iy 1 (4-30) 

+ 



M(cosh k - 1) - iQ 2 M(cosh k - 1) + iQ 2 

where Qi = Msinh/c, (fi(0) is the form factor function (|2.36|) and the constants mi 
are given in (|2.38|) . Since for small k 

(j){iTX + k) = --, (4.31) 

the 1-particle contribution, for fixed Q 2 as function of Q2, has poles at Q2 = iy^Q 2 
with residue —^5n. 

Assuming that no other singularities are generated by the higher intermediate 
states away from the imaginary Q2 axis, from the usual Cauchy integral we conclude 
that for a circular contour centered at the origin with radius y/Q 2 < R < Q 2 /(2M) 

1 /dCT7,(Q 2 ,C) . N ( ~ .2M 1 (AVJ . 

- ■ r kN (Q ) - t —5 llT7 =—— , (4.32) 



2-KiJ C N+l ' V ' Q 2 (^) N+1 

where ti-n{Q 2 ) are the coefficients of the Taylor expansion 

00 

ri(Q 2 , Q2) = t I;n(Q 2 ) (iQif ■ (4.33) 

N=0 



12 



Now the structure function moments can be computed in the usual way by calcu- 
lating the Cauchy integral along the deformed contour around the cuts. In this way 
we obtain expressions for the moments: 

M iiK (-Q 2 ) = M ^ (<fl + f fc^-^j} • (4.34) 

The spin function f has no 1-particle contribution, and so assuming no further 
singularities apart from the cuts we obtain for the moments 

~ / O 2 \ N+1 

M l . N (-Q 2 ) = M\^j f l:N (Q 2 ), (4-35) 

where 

oo 

n(Q 2 , Q 2 ) = £ n, N (Q 2 ) (iQ2) N ■ (4.36) 

N=0 

Note in the equations above for / = 0, 2 N is even, positive and for I = 1 N is 
odd. It remains to extract information on the Taylor coefficients t/ ; 7v,t/ ; at from the 
operator product expansions. 

4.3 Operator product expansion for the spin field 

Starting with the spin field, the connected part of the time ordered product can be 
expanded as: 



T E (f cE (yi,y 2 )1> d (0,0)) 



c 

oo 



E^ )C ^V) + E E ^'V) {^ )C V +^ )C V} 
j=i 



(4.37) 



where y± = =Fyi — ii/2- Employing a basis of hermitian operators 

A^) cd ^ = A$ cd , qW)^ = gW) cd (4.38) 

and using Poincare symmetry, parity and CPT invariance we have 

^{J,l)cd = v B{ j il)cd v ^ A{ i )cd = v A{ i )cd y ^ (4 ^ 3g) 

where V is the parity operator and 

iF(y 2 ) = tSV), i { J' l) \y 2 ) = tW)- (4.40) 
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Further we define the matrix elements B^ ' as in appendix [E] 

./ 



W) 

(a, 6\ B { J' l)cd \b, 9) = ^yej " Pf> cd (4.41) 
and we find 

flWO* = flWO = (-1) J + l flWO . (4.42) 
The "twist" of the operator is defined as 

4 J, ° = dim (b^ *) - J (4.43) 



and the minimal possible twist value is zero. The contribution of these operators 
dominate for large momenta and we have 

M/;7v(-Q 2 )«T^r)W)(Q 2 ), (4.44) 

where 

f, W \Q 2 ) = (Q 2 Y +1 (^) J /dV Q V J 'V), (4-45) 

n w \v 2 ) = E^7W)- (4-46) 



4.4 Operator product expansion for the current 

Using hermicity, Poincare, O(n), parity and CPT symmetries and current conserva- 
tion we can write 



i=0,2 

(4.47) 



where the dots indicate that we have omitted total derivative operators since they 
would not contribute to the diagonal expectation values. Otherwise the set of op- 
erators appearing here is as in Q4.37JI and the coefficient functions H^l-^ take the 
form 

hII-Av) = - y -y(y + (y 2 ) , (4.48) 

H^(y) = H%Jy) = (y + 1) 2 V« (y 2 ) , (4.49) 
H { l^y) = -^y(y + 1)V® (y 2 ) , (4.50) 

y_ 
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where V® (y 2 ) (I = 0, 2) are real functions unique up to co ^ st ' , and 

,_d 

dy 



y = y 2 4-2- (4-51) 



Similarly 



= —Sji - y-y J +-\y + J)(y + J + i)vJ J>,) {y 2 ) , (4.52) 
y+ 

Hi%(y) = H { _%(y) = y{(y + l)(y + J + 1)F «<) (y 2 ) , (4.53) 

= - v -^y{y + i)vS J ' l) {y 2 ) , (4.54) 

where V^ J '^ (y 2 ) (J = £) are real functions unique up to co ^ st ' 5 j\ and c w are real 
constants. Finally 

H% 1 1(V+,V-) = H$l{y-,y+), (4.55) 

where Jx = —fx for the light-cone index fx = ±. 

As a consequence of the asymptotic freedom of the O(n) model for small y 2 

vS l) (y 2 ) ~ \yf~ 2 , (4.56) 
(y 2 ) ~ \yfi' l) -\ (4.57) 

where tH\t^'^ are the twist of the corresponding operators. 

4-4- 1 Fourier transformation 
Introducing 

xU(y)=y 2 vS l \y 2 ), (4-58) 

xW\y)=y 2 y{vW\y 2 ), (4.59) 

X { f\y)=y 2 y J _vS J ' l Hy 2 ) (4.60) 

the current operator product in Fourier space can be written 

d 2 Qe^(a, 0|T E (j^ E (y) J e J (0)) % 0) = - £ Bf^ f E^{Q)X^{Q)A^ 

1=0.2 
U) 

Y^Rf^^E^Q) [xW) (Q)+ fi J * (Q)] +mc a;( 5 J1 ^(Q)} (zMfj^W) 

(4.61) 
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where the reduced matrix elements A$ , B^f' 1 ^ are real and E fMl/ (Q) is the transversal 



tensor 



E (U ,(Q) = QvQu-Q%v. (4.62) 

The complete expression (|4.61|) . although conserved in coordinate space, is not 
transversal because of the anomalous terms proportional to the constants c w . These 
are multiplied by the tensor K^ U (Q) with components 

K ++ = -^-, K - = 7^-> K + _ = K_ + = 0. (4.63) 

It is not quite trivial to see, but easy to check that 



k,AQ) = - E ^Q)^i + ^7 - ^ . Cm," = ±) • (4-64) 



Thus is transversal up to the last two terms, but these correspond to con- 
tact terms in coordinate space. Dropping these "seagulls", the coefficient of the 
transversal part in Fourier space becomes 



n{Q) = -\^X${Q)AV (4.65) 

QlQ 2 



5E (=^) V> {n j m+t J ' l) (Q)-^ Sjl ^- c 



This can alternatively be written as 

n(Q) = -\Y.^ ) (Q) A( £ > -\Y.{ [(2MQ_) J + (2MQ + ) J ] (4.66) 

)V™^ )B ^} + M,^*_ a , lE c„ B <»., 



AQ 2 

where 

W«V) = y VWV). (4.67) 

Note that the anomalous contribution to (|4.66j) is regular on the imaginary Q2 axis 
hence does not contribute to the structure functions. Let us also define 

e w V) = £^V)^. ( 4 - 68 ) 
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If we now compute 

J 



£W\Q 2 ) = £ (Q 2 ) J+1 ( JL) W™(C?)bM (4.69) 

we see (using asymptotic freedom) that the coefficient functions behave as (Q 2 ) " , 
up to logarithmic corrections. We will keep the contributions of the leading (twist 
0) operators only. Note that for J = I = 1 the only twist operator is fi^' 1 ^ ab = J<^> 
with 

#J 1 ' 1) =4. (4.70) 
From ()4.66|) we obtain for the Taylor coefficients 

^ (< n_ W) _*££i^_, (4 . 71) 



where (up to higher twist contributions) 



3 



-iMQ 2 )^-\¥^i m) (Q 2 ) ■ (4-72) 



2 (Q2)^v+i 

Inserting this in Q4.34|) we obtain: 

!2 \ W+l 



1 



M, ;N (-0 ) = M (i^j |^ (Q 2 ) + (1 - 2. ci)j 

(4.73) 

Later we will see that c\ = l/2n. Thus these two subtle effects cancel each other, 
and so the final formula coincides with the naive one: 



Q 2 ) = (Q 2 ) ■ (4.74) 



Up to now we have related the moments to the Taylor coefficients which we see 
are determined by the structure of the OPE. But to get quantitative results at this 
stage we need more dynamical input. This can be supplied by analyzing the OPE 
in the framework of renormalized PT, which is the topic of the next section. 



Note for large Q the anomalous term dominates over ti-.n- 
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5 Perturbation theory and operator product expansion 

We consider the 0(n) sigma model Lagrangian 4 

1 n 

and work in D = 2 — e dimensions using dimensional regularization. Renormalized 
fields and coupling g are given by 



2 

S a = z l/2 j Z = 1 - + • • 

2 e 2 ^ 7 1 2f3 g 2 

9o = V 9 z i , Zi = l h 



(5.2) 



We denote the usual renormalization group (RG) derivative by 



where the dimensional regularization beta function is 



P(jg) = -^g + P(jg), P(g) = -Pog 3 -Pi9 5 -... (5.4) 



and 



The RG A-parameter in the MS scheme satisfies ^Aj^g = and is written 

A Ms = ^ e/(9) > ( 5 - 6 ) 

where 

with 7 = ln47T + r'(l); (note ^ = ^). 

The spin field $ a differs from the renormalized O(n) field S 1 ^ only by a finite 
renormalization: 

$ aE = n n S% (5.8) 



4 In practical computations in infinite volume one usually adds a coupling to an external field 
(S n - 1) to serve as an intermediate IR regulator. For IR finite quantities the renormalized 

9(1 

external field ha. = ho\/~Z/Z\ is set to zero at the end of the computation 
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and solving the RG equations for the vacuum two-point function the standard way 
we find 



n-l 

27T \ "- 2 



where the constant C n is that appearing in Eq. (|4.12|) and p(g) is the solution of 

where 7(g) is the anomalous dimension of the spin field: 

7(5) =V\nZ = lQ g 2 + ■■■ , 1q = VL—_ (5.11) 
The integration constant in (j5.1U|) is fixed by requiring 

{<?)**> {l + 0(g 2 )} . (5.12) 



e - P (g) 



5.1 Zero twist operators 

We now introduce a basis for zero twist operators composed of an even number of 
spin fields. For isospin I = we write 

K (n umi )-(n k , mk ) = ±_ (Qniga! . #£15*1) . . . (#*5*»A . d^ k S ak ) , (5.13) 
% 

where we introduced the notation 

d± = ~(id 2 Td 1 ) = ~(d oT d 1 ). (5.14) 
It is very important to notice that a complete basis can be chosen such that 

nj,mj>l, i = l,...,k, (5.15) 

which can be achieved by using the identity 

m-l / v 

S a Qrnga = _ V- m gm-i^a . gi^ga _ ( 51 g) 
2 <=1 V*/ 



The spin of the above operators is Yli=i ( n « + m «) = ^ ; whereas the mass dimension 
is J — ke, i.e. the operators are of zero twist only in exactly two dimensions. For 
/ = 1, 2 we can define the operators 

j^ab(n ,m )(n 1 ,mi)---(n k ,m k ) _ quo go, _ gmo j^(m,mi)—(n k ,m k ) (5 ^\ 

% 
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with spin Y2i=o ( n * + mi ) = ^ > an< ^ dimension J — {k + l)e. For ! = 1 we have to 
antisymmetrize the a, 6 indices, whereas for I = 2 we take the symmetric, traceless 
part. Now 

< no < rriQ , rii, mi > 1 i = 1, . . . , k (5.18) 

and correspondingly there are three types of operators: 

type I: 1 = and I = 1, 2, no > . 
type II: I = 1, 2, no = 0, mo > 1 . 
type III: / = 2, no = mo = . 
It is now straightforward to calculate the potentially divergent matrix elements 

of these operators. At one-loop order, after wave function, charge and mass renor- 

malization we find 

1 = 0: all matrix elements finite 
( 9 2 \ 

1 = 1: type I : 1 x lowest order 

type II : only type I operator matrix elements 
( 9 2 \ 

1 = 2: type I : 1 x lowest order 

V 27re / 

( 9 2 \ 

type II : 1 x lowest order + type I 

v ™j 

( 9 2 \ 

type III : 1 x lowest order + type I, II 

V Trey 

5.1.1 Operator product expansion at tree level 

The leading terms of the OPE in perturbation theory are simply given by Taylor 
expansion: 

T E (\s a (y)S\0) 



9o 



c 

oo 



, 1-1 (5 - 19) 



yi 



up to higher twist operators. The operators appearing in the sum over J can be 
written as a sum over operators of definite isospin: 

— d J + S a -S h = ^O ( ( o)° a \ (5-20) 
ff o i =0 

°( ( o J )° ab = \p? b;cd diS a ■ S b . (5.21) 
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Important operators with isospin 2 are 



T (°) = ^ ( SaS " ~^ Sab) - (5 ' 22) 



For isospin 1 we have the currents 

J So) = \ (s a d,S b - S b d,S a ) , (5.23) 

and for isospin we have the energy-momentum tensor 

T M0) = ^ (d»S a d v S a - 1 V drFdvSr) . (5.24) 

In terms of these, the leading operators of the OPE can be written as 

(j,2)afe = + type j operators j ( 5> 25) 

(ji)afe = Igj-ijal ' + type I operators , (5.26) 



(0) 2 + +(0) 

,(2,0)a 6 _ 

n 



0(0) = -- T ++(°)- ( 5 - 27 ) 



5. J. 2 Renormalization of the zero twist operators 

We will now denote by B^^ ab the zero twist operators introduced in the preceding 
section. Here a is a multi-index: it includes the operator type I, II, III and possible 
further indices. In D space-time dimensions the mass dimension of B^^ ab is 

J — e cfa' 1 ^ . The corresponding renormalized (finite) operators of mass dimension 
J are: 

B W)a b = £ Z W) ^ B W)a b j (5>28) 





where the operator renormalization constant matrix is 

t 
e 

We now distinguish the types of operators by writing their multi-indices 



Z ( a f = <W " 9 T ™lf + ■■■■ (5-29) 



B 



(J,l)cd 
a(0) 



yj(J,l)cd 

„(J,T)cd 



for type I 

for type II (Z = 1, 2) (5.30) 
for type III 
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In this notation the one-loop results of the previous subsection are 

(J,l) _ 1(1 ~ 3) W W)_ W W)_ Q (531) 

w ii = ~\ 5ab - (s - 33) 

Using the inverse matrix W: 

E z if < l) = ^ > <f = *<# + t w % ] + ■ ■ ■ ( 5 - 34 ) 

we define the anomalous dimension matrix 

*° = E (*™&'° " , (5-35) 



which is finite and is given by 



uW) = a 2 W W> 
pa y w 



{d W)_ d W)-l) + {g A ) (5.36) 



pa H w pcr \ CT ^p 

to leading order. 

We would like to go to a basis where the leading anomalous dimension matrix 
is diagonal. This basis is easily found due to the triangular structure of the leading 
anomalous dimension matrix. The renormalized operators in this new basis are 
denoted Bu?'^ cd , where u> = a,A and B^ ,,l ^ cd are (as before) the renormalized type I 
operators and Z3^'^ cd are operators of type II and (for I = 2) of type III mixed with 
lower type operators. In this basis we have 

«#° = ^^W + OG, 4 ), (5.37) 
S A B9 2 + 0(g 4 ), (1 = 1,2), (5.38) 
0(g i ), (1 = 1,2). (5.39) 

r ab 
r (0) 



(J,i) Kg ~ Q x 2 , u U 

W) 1 ~ 1 X 2 

".»/; = — ' 

V aA - v Aa ~ ^xa ) i V — *-■>")■ 

We also note that the canonically dimensionless I = 2 operator is multi- 



plicatively renormalized: 

r^=Y^T ah , (5.40) 
since there is no other operator with the same quantum numbers to mix with. Here 

Q 2 

Y = !- — + ■■■ (5.41) 
ire 
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leading to 



a 2 

V InY = ^ + 0{g 4 ) . (5.42) 

7T 



It is clear that for the / = operators we can chose the diagonal basis so that 

o (J,0)cd = B (Jfi)cd (5>43) 

and for / = 1, 2 the A = 1 operators so that 

B W)* = d J T c* =fl -eYB[ J > 2)cd , where B[ J ' 2)cd = d J + T cd , (5.44) 
B[ J ^ cd = dl' 1 J c + d {0) =^B{ JA)cd , where B[ J ' 1)cd = d^ 1 J? . (5.45) 
We then have 

n (J,2)cd _ „(J,2)cd ST^ X (J,2) n (J,2)cd 

^(0) ~ D l(0) + Aa a(0) ' {0-40) 

a 

a 

Finally for (J, 0) = (2, 0) there is just one operator and we have 

B (2,0)cd = _^ T++(()) = ^ B (2,0)c d ^ where B (2,0)c d = J_1 T++ ^ (5 . 48) 



and therefore, the 1-particle matrix element is known exactly: 

B[ 2fi) = 2 . (5.49) 

5.1.3 The operator product expansion in perturbation theory 
In bare perturbation theory we have 

T E \(s a (y)S b (0)) = 

E4!»" 2 > + e e *it»y> k$°v + s s>°v} . <5,50) 

Z,u; J=l £,u; 

which, after renormalization, becomes 

T E (s a R (y)s b R (o)) c = 

OO (EC r-|\ 

E ^ )ab (y 2 ) + E E ^ J, ° (^ 2 ) K* /)afc w+ + ®« l)ab y- } > 
i,u j=i i,u 
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where 

*W) = S ^ E k m(y 2 ) ^ ww> (5.52) 



Z ^ p(o) 
p 

which satisfies the renormalization group equation (RGE) 

(P + 7 ( 5 ))A^)=E^ J '°^ - (5-53) 

Finally the original coefficient functions 7^ of ()4.37j) are related to the renormal- 
ized coefficients k^f' 1 ' by 

4 J ' l) = nlkW\ (5.54) 

and so ()4.46|> can now be written as 

^V) = ^E^'V)^ J ' - (5-55) 

U) 

The perturbative expansion of the renormalized coefficient functions is 

kW\v\y\,g)=g 2 K^+gUi J ' l \»\y\) + 0(g 6 ). (5.56) 
We already computed the leading (tree-level) terms: 

K^ 0) = ±8 ul , (5.57) 

4 J ' 1} = ~ ^i, K[^ = ^\ (5.58) 

Kf 2) = j { 5 A1 , KW = ±\W. (5.59) 

A building block used in the solution of the RGE is the matrix Uua- l \g), which 
is a solution of the matrix differential equation 

Pig) ^ W/Hg) = -E "#°G0 ^GO- ( 5 - 6 °) 

9 p 
If we have such a solution and its matrix inverse U satisfying 

Z ^uW\g)UW\g)=5^, (5.61) 
p 

we can build the RG invariant coefficient 

G W) = #0) kW\l,g) uffHg) (5.62) 
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and the RG invariant numbers 

V<^=J2u(f/\ 9 )BW). (5.63) 
p 

In H5.62JI the running coupling g is defined as the solution of 

f(g) = f(g) + Hfi\y\), (5.64) 
which, for small \y\, has the asymptotic expansion 

2p f = A + cA 2 + 0(A 3 ) , c= i(r'(l)-ln7r), (5.65) 

where the effective coupling A is defined by 

1 1 ~ 2e r 'W 

■=r + In A = In- — . (5.66) 

A n-2 Aj^lyl 

Putting the building blocks together we have 

n-l 

*W) = (^Y~ 2 nC^G^vW). (5.67) 



Note that if with some Y 



B[ J $ cd = »- e YB i 1 J ' l)cd , (5.68) 



which is the case for 1 = 1,2 and also for (J, I) = (2, 0), then 

v[ J f =VlnY6 lu; . (5.69) 

5.1.4 Solution of the matrix problem 

In this subsection we will omit the upper index W> and use matrix notation. We 
want to solve 

0(g)^-U(g) = -u(g)U(g), (5.70) 

which is (|5.6Uf) in this notation. We know that in our basis 

u(g) = 2(3 Ag 2 + O(g 4 ), (5.71) 
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where A is a diagonal matrix with diagonal elements 
1 = 0: A a = 0, 

1 = 1: A A = 0, Aa = ~ 2' ( 5 - 72 ) 

2 . 1 

Z = 2 : A A = -, A a 



n-T n-2 
Using the expansion 

^4 = _^_ 2 y ^AW, (5.73) 

we can take the Ansatz 

f>( 5 ) = [l + J R(<7)] 5 2A (5.74) 

with 

oo 

%)=^ff 2s K (s) , (5.75) 

s=l 

and put it into (|5.7U|) . We get 

s-1 

+ [i? (s) ,A] = +^ A^flW, a = 1,2,..., (5.76) 

p=i 

which has a unique recursive solution unless A w — A, T = s occurs for some w, a and 
s. In our case this is possible only for s = 1 and only if n = 3. For n = 3 we thus 
take the modified Ansatz 

= [1 + R(g) + In 5 2 R(g)} g 2A (5.77) 

with 



s=l 

In this case we start with 







oo 

-E 


5 2s J2« 


r2 = o, 




J 


s(i) 

-"-AB 


_ 4(1) 


o(l) _ 1 


,(i) 


p(l) 



(5.78) 



' = 1: fl» = 4u(. *£=0, «S = 0. ^ = 0, (5-79) 

?(i)_n dW- 1 ^!) D (i) f ( 5 _ 80 ) 
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and 

' = 2: R { H = A^l R%=0, R$=0, flW =0> (5 . 81) 

p(l) _ n tfW - 4W _ 4 (1) o(l) _ 1 4 (1) /t:oo\ 

and after that there is a unique, recursive solution of the system 

s-l 

sR (s) + [R {s) ,A] = ^ A^-d , (5.83) 
p=i 

s-l 

s i? (s) + , A] = A« - flW + ^ A (S ~ P) R (P) (5-84) 

P =i 

for s = 2, 3, . . . 

Note that from the recursion relations it follows that 

fiftl = for 1 = 1, (5.85) 
R { *\ = for 1 = 2. (5.86) 

We also note that because v\ w is proportional to 8\u 

^t )=0 and R^~ S iu> (5-87) 

and therefore 

Uiu = — 5iw U lw = w5i w , (5.88) 

w 

where w is the solution of 

0(gW(g) = y(g)w(g) , (5.89) 

where y(g) = T> In Y (g), the coefficient occuring in Eq. 1)5.69 J) . This has the following 
important consequences. 

Vi = for I = 2 , and for I = 1, J > 1 , (5.90) 

Vi = Si for Z = 1, J = 1 , and for I = 0, J = 2 . (5.91) 

5.1.5 Leading terms in coordinate space 

Using the results of the preceding subsections we can calculate the leading terms in 
the short distance expansion of the functions (]4.46[) . We find 

ri U0) = jy/ (J) A-^{l + O(A)} , (5.92) 
^(1,1) = f[V\-^2 jl + 0(A)} , (5.93) 
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where 

^ = (5.94) 
A (1) = -^BW. (5.95) 

For the case I = 1, J > 1 we have to distinguish between the cases n > 3 and 
n = 3. In the former case 



= const. + — f[ J) {l + O(A^)} , (5.96) 



where 



and 



2 

In the n = 3 case we have 



// = const. + In A 1 1 + 0(A)} , (5.99) 



where 

I 

7T 



a (j) = - E 4 J) ^ , 4 J) = E ^ *21 W1) • ( 5 - 10 °) 

Finally for l = 2we find 

W2) 



where 



and 



// = const. + — f { 2 J) A {l + 0(A^)} , (5.101) 



2n-3 

ti J) =(-^) n ~ 2 nC^^VyV, (5.102) 



4 .) = ^X-W + £ A W2) ^Wa) + j< ~iJ» (5 . 103) 



Using the coordinate space results above and the asymptotic formulae of Ap- 
pendix[GJ we are now in a position to derive the results on leading large momentum 
behavior of the spin structure function moments given in subsection 4.1. 
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5.2 OPE for the currents 

It is straightforward to calculate the leading operator product coefficients in pertur- 
bation theory: 



^SV)=-5, (^T^ i+0(1) ' J - 2 (5 - 104) 

W S^ 2 ) = 2^(^2)^^ fe + 0(1) ' J " 3 (5 - 105) 

ywl$(y 2 ) = ~ (l " |) + O ( 5o 4 ) ■ (5-106) 
where the operator associated with (|5.1(J5j) is 

B k(0) ab = ^2 ( d +~ lsa ■ d + sb ~ di^S* ■ d+S a ^j . (5.107) 

Eqns. (|5.1U4|) and ()5.1U5|) can be obtained by tree-level perturbation theory, while 
the results necessary to write down the one-loop formula (j5.106[) can be found in 
jllj . Also the results of |3lllj show that the coefficient in (|4.71|) is given by 

ci = ^. (5.108) 
Using renormalization group improved perturbation theory we can write 

e wo (y2) = ^ r wo (y2)y w) ) (5 . 109) 

where, as in (|5.(i3j) . 

yW = £ff«O(0)flttO (5.110) 
p 

are renormalization group invariant constants and 

rW\y^ = Y,wW\y^uW\g) (5.1H) 

p 

are renormalization group invariant coefficient functions. Putting everything to- 
gether we arrive at the results already given in subsect. 4.1. 

6 Structure functions for n = 3 

In this section we consider the case n = 3, where it is possible to compute the 
structure functions accurately in the whole range of x for a given q 2 . The case n = 3 



1 
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is rather special for various reasons e.g. the spin and current 2-point functions 
exhibit in this case very similar properties and there are miraculous scaling relations 
|H] which relate them 5 . In the S-matrix bootstrap approach its distinguishing 
feature is that it is the model for which the r-particle form factors can most easily 
be obtained explicitly. They take the simple form 



n, ■ 



where 



9, 



-7T 



3r/2-l 



n 



l<i<j<r 



„ -tanh 2 — 
(9(2vrf -9) 2 



(6.1) 



(6.2) 



(6.3) 



and the reduced form factors g% br are polynomials in the rapidities. There are 
well defined recursive procedures for computing the form factors, the only practical 
limitation being that they become very involved. So far the record we have achieved 
is the 7-particle form factor of the spin field ^Hj > already its algebraic expression in 
MAPLE involves many megabytes of storage. Fortunately for the structure functions 
we only require sums over bilinear factors of the form factors which are computa- 
tionally more manageable. In correspondence to (|6.1|) we define reduced form factor 

(r) 

squares j\ through 



\y r+1 (iTrJ 1 ,...J r )\ 2 tf ) 



1 



-7T 



3r+l 



n B(u jk ) 

l<j<k<r 



(6.4) 

jl r \p u ...,P r )(6.5) 



where we have introduced two new functions 

A(0) = ^{m-ef- 



B{6) 



{6 2 + vr2)2 tanh 4 | ' 
02(02 + 4 ^) tanh 2' 



(6.6) 



(6.7) 



The reduced form factor squares for r = 2, 3, 4 are given in Appendix EI F° r r > 4 
the expressions are too lengthy to exhibit in print; the results for r = 5, 6 can be 
obtained in the form of files from the authors. 

For r = 1 we then have (noting that sinh = k for x = 1), 



w 



(l) _ rniTT 



A(b)6(l-x) 



(6i 



where the factors m; are given in (|2.38|) . 



see also the OPE in sect. 4 
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6.1 Case r = 2 

For the case r = 2 the delta-function constraint in the integral representation is 
simply solved and we obtain the analytic expression 

6 2 

wf\q\ x) = 9{u - 2)—^==A{h + 0/2)A(A - 0/2)Bfa)C, (A, 0) , (6.9) 
where the kinematic variables u,A,<f> are given by 

2 cosh | = w, (6.10) 

w = W/M, T^ 2 = (p + g) 2 , (6.11) 
coshA = (M 2 + pg) / (MW) , (6.12) 

and 

C*(A, 0) = jf } (A + 0/2, A - 0/2) . (6.13) 
Using the expression for in (|H.2|) - ()H.4|) we have, 

C = 8vr 2 + 4A 2 + 30 2 , (6.14) 

d = -2vr 2 + 2A 2 + ^cf> 2 , (6.15) 

3 

C 2 = 2vr 2 - 2A 2 + -<f) 2 . (6.16) 

Despite its relative simplicity, this case exhibits many features in common with 
higher r. The structure function approaches its asymptotic values very slowly e.g. 
for q 2 fixed 

wP'(o 2 ,x)~e; rr t r- for x — > , (6.17) 

with eo = 1, ei = e2 = 1/4 (consistent with the small x behavior derived in sect. 3), 
while for — q 2 — > oo , x fixed we have 

(2) 7r 6 x,4(-ln(l-x)) 

~ e/ 8(l-x)(ln(-,VM 2 )) 2 - (6 ' 18) 

6.2 Results for the entire x range 

Just as for the 2-point functions [2] we find that for a fixed — q 2 only states with 
a limited number of particles contribute significantly. To appreciate this better we 
consider the sum of the field and current structure functions, which is a rather 
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Figure 1: Approximations to x(wo + wo) as functions of < x < 1 for —q 2 /M 2 — 100. Curves 
correspond to sums up to and including 2, 3, 4, 5, 6-particle intermediate states. The last 3 curves 
are indistinguishable on this scale. 

peculiar thing to do in general, but which is in fact rather natural in the special 
case n = 3. Figs. 1 and 2 illustrate how the structure function x(wq + wq) is built 
up from states with increasing particle number for the cases —q 2 /M 2 = 10 2 and 
—q 2 /M 2 = 10 4 respectively. We see that the higher states contribute very little 
and that we obtain nearly exact values for the structure functions for all values 
of —q 2 /M 2 < 10 5 by including only intermediate states with < 5 particles for the 
current and < 6 particles for the spin field. 

In Fig. 3 we plot xwo(q 2 , X^j clS 9j function of log 10 (— q 2 /M 2 ), for a selection of 
x-values 6 . The function increases as — q 2 increases for all values of x in this range 
and seems consistent with Bjorken scaling as mentioned in sect. 4. 

6.3 Threshold behavior 

Note that in Fig. 2 we have cut off the plot at x = 0.95. This is because near x = 1 
the function develops a big bump with a peak ~ 70 which, if included in the same 

6 For this model we prefer to show this rather than the typical HERA plot where one adds 
— log 10 (x) to separate the a;-values, because the latter would obscure the —q 2 variation which is 
rather small compared to the variation of — log 10 (a;). 
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0.2 0.4 0.6 0.8 

Figure 2: Approximations to x(wo + wo) as functions of < x < 0.95 for —q 2 /M 2 — 10 4 . Curves 
correspond to sums up to and including 2,3,4,5,6-particle states. 

plot, would obscure the features we wanted to show there. The behavior of the 
a-model structure functions near x = 1 is indeed rather involved. For a fixed — q 2 
the contribution to the structure function from the r-particle state vanishes for 
x greater than some threshold value 

x r (-q 2 ) = [1 - (r 2 - 1)M 2 /9 2 ] _1 • (6.19) 

The big bump in x(wo + wq) referred to above is at this value of —q 2 /M 2 = 10 4 
practically entirely due to the 2-particle contribution. For this contribution: 

wf^ ~ Ei{q 2 ) \J xi{— q 2 ) — x , x —> X2 , — q 2 fixed, (6.20) 
wf ] ~ Fi{x)/\n 2 {-q 2 /M 2 ), -q 2 -> oo , x fixed, (6.21) 

where E[,Fi are some (known) functions. The bump arises because Fi is quite 
singular near threshold, Fi ~ [(1 — x) ln 2 (l — a;)] -1 . The analytic behavior as x —> X2 
sets in only extremely close to threshold e.g for —q 2 /M 2 = 10 4 the position of 
the peak of the bump is at x = 0.99954 whereas the function vanishes at X2 = 
0.99970. At —q 2 /M 2 = 10 4 the 3-particle contribution also has a bump but it is less 
pronounced; (peak value ~ 2.5 at x ~ 0.9953). We conjecture that the threshold 
behavior of in the 0(3) model is (x r — x)^ 2 ^ 3 ^ 2 . 
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log_10(-q A 2/M A 2) 

Figure 3: xwo(q 2 ,x) for various values of x — 10 -1 ^ 5 . 
6.4 Moments 

For r = 1 the moments (j2.29|) are simply given by 

4 

M^(q 2 ) = ^-A(2 asurtus) . (6.22) 

The moments with N > 1 are quite simple to evaluate numerically, but for N = 1 
some care must be taken to obtain accurate results. 

The problem arises already for r = 2 where we have integration just over u\. 
The ^-factor together with the x 2 factor in the integrand is 

Ff ] = A{h)A{h)B{u{)x 2 . (6.23) 

Now for u\ very large ~ (~Q 2 + M 2 )e~ Ul is exponentially small and so also 
x ~ — q 2 e~ Ul . Noting i) for u\ large (3\ ~ u\ and ii) A{9) ~ 16/(tt 4 9 2 ) for 6 — > we 
have for large u\ 
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2 4 6 8 10 

log_10(-q A 2/M A 2) 

(r) 

Figure 4: Contributions Mq. 2 for n = 3 from r = 1, . . . , 6-particle states. The upper full line is 
their sum. The dashed line is the perturbative expansion of Mo-2 + Mo-2 = 1 + A up to and including 
terms of order X(q 2 ). 



where 

A B (u) = A(u)B{u) = ] , (6.25) 

The integral over large u\ gives a sizeable contribution because the integrand decays 
only as u^ 2 . The integral is thus broken up into two parts where for the large 
u\ region the substitution (|6.24j) is made and there computation of exponential 
functions of large argument are not necessary. 

For the case of higher r the procedure is similar. Here the integrations over 
ui, . . . ,u r _2 can be done safely by introducing for them (large) cutoffs (and moni- 
toring the dependence on them), since the integrands are exponentially suppressed. 
But for the integration over large u r _i the integrand is not exponentially suppressed 
and in this region one replaces the corresponding ^-factor by 
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2 4 6 8 10 



Figure 5: As for Fig. 4 but for the moment I = 1, N = 1; here the PT result is 1 + 0(A 2 ). 



In Figs. 4 and 5 we plot the separate r-particle contributions Mq. 2 and M}.-/ 
respectively; some corresponding numbers are given in tables 111 Appendix [IJ 
They are typically bell-shaped (except for r = 1) and perhaps obey scaling relations 
similar to those of the spectral functions examined in ref. 9 . The figures show how 
they build up the sum of moments Mq ; 2 + Mq-2 and M\-x + M\.\. Using the exact 
ratio of the mass to the A-parameter Eq. (|4.18|) . we also exhibit the perturbative 
results up to and including terms of order \{q 2 ). The agreement of the summed 
terms and PT is impressive for —q 2 /M 2 ~ 10 5 . For values of —q 2 /M 2 >~ 10 6 
contributions from states with > 7 particles must be taken into account. Note we 
have also included the contribution of the one particle states in the sums; these 
tend to improve the agreement at lower values of — q 2 and fall asymptotically as 
Mg} ~ mi7T 4 /[Aln 2 (-q 2 /M 2 )]. 
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7 Sigma model structure functions in the 1/n approximation 



7.1 The spin field structure functions 

In the framework of the 1/n approximation the spin amplitude (|2.19j) has an expan- 
sion of the form 



x 1 



Y 1 ab ^ d (p,q) = Y J ^ T . S\ ] 5 ac 5 bd + S [ z ] 5 ad 8 bc + Sl ] 5 ab 5 cd , (7.1) 

r=l n 

and so for the various isospin channels 

oo 1 

m(q 2 ,x)=J2-4\<l 2 ,x): (7-2) 

r=l 

with 

fijH = fif ] + 4 r] + 4 r+1] , (7.3) 

w [ ; ] = s [ ; ] -s?, (7.4) 

= + (7.5) 

The Feynman rules for the 1/n expansion of the <r-model has been described in 
many places (see e.g. ref. |14| and references therein) and will not be repeated here. 
We only mention that the diagrams involve the bare propagator of the fundamental 
spin field, and the bare propagator B^ 1 of an auxiliary isospin scalar composite field, 
which we call A, which is the inverse of the scalar 1-loop integral given in App.ljJ 

In leading order 1/n the only contribution to the scalar structure function is the 
tree diagram involving A exchange in the "s-channel" ; one thus gets an amplitude 
proportional to the imaginary part of B : 

5™ = 4 1] =0, (7.6) 
oW a at o\ -q 2 M 2 sh(j) 

where <fi is as in ()6.1U|) . Note we already anticipated = by starting the sum 
over r at 1 in (|7.2|) . 

In the limit of small x we have 

^i 11 ~ T\ aJT/ 2\2 ' 7~2 > x ~^ ' -I 2 fixed ' ( 7 - 8 ) 

(1 - M 2 /q 2 ) 2 xln z x 

consistent with the general result ()3.15j) for / = 1,2, the scalar Adler function in the 
leading order 1/n being just Aq(z) = (l-\-M 2 /z)~ 2 . Note the limit is approached ex- 
tremely slowly. E.g. denoting the asymptotic function on the rhs of (|7.8|) by .Sasympt 
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one has sj 1] /5 asympt = 0.245,0.441,0.638,0.907 for x = 10" 5 , 10~ 10 , 10~ 20 , 10 -100 
respectively. 

One also observes the limit 

_rii 2n(l — x) 1 9 _ . . 

S\ 1 i L . — - q 2 _> oo x fixed , (7.9 

1 x \n 2 {-q 2 /M 2 ) K ' 

and the threshold behavior: 

~ q i^2\9 — V^2 - x , x -> x 2 , -<? 2 fixed . (7.10) 
7r(l — q z /M Z ) Z X2 

We caution that the limits n — > oo and x — > x 2 may not commute. 
The moments (j2,28jl have an 1/n expansion of the form 



oo 1 

%tf) = £^M{ r i(? 2 ), (7.11) 



n 

r=l 



M /M^ 2 ) ~ 777^7— 7777^7 for ^ > 1 , 2 = 1, 2, (7.13) 



One then shows (e.g. numerically) that for — g — * oo 

N(N -l)\ii 2 {-q 2 /M 2 ) 

consistent with the results (j4.15|) and (j4.16j) . The l/ln(— q 2 ) behavior for the N = 1 
moment comes from the singular behavior at x = 0. 

So far we have only obtained the leading order for the isospin I = 1, 2 channels. 
This is because Eq. (|7.3|) shows that to obtain the leading order 1/n approximation 
for I = 0, one has to take into account also the amplitude . To this amplitude the 
only contribution comes from a diagram with two A-exchanges in the "t-channel" 



4 21 = 27r( _ g2 4 f M2)2 Im / H^fc + ? E )^ - P E )B(k)- 2 , (7.14) 



where is the Euclidean bare spin propagator 

d(k) = (k 2 + M 2 )~ x , (7.15) 

and q E = (iqo,qi) and similarly p E with p E2 = —p 2 = —M 2 . Using the spectral 
representation of B{k)~ 1 and the cutting rules 7 (see AppendixQJ we get 



7 One can initially introduce a UV cutoff in the spectral integral and remove it after invoking the 
cutting rules 
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where 



k 2 



KM* 



± 



x 2 (-q 2 ) 



V 



R 



\l-y)(l + 2x 2 y 2 /R) 
l-y + 2y/R 



2M 2 x 2 (-q 2 ) ' 



Again the small x limit 



-[2] 



2tt 



1 



, x —* , — q 2 fixed . 



' (l-M 2 /q 2 ) 2 xln 2 x 
is as expected. For the large q 2 and threshold behaviors we have 
,[2] 2n 1 



8$ - 



^ 3 ~ 



x(l - x) ln 2 (-g 2 /M 2 ) ' 

8ir(-q 2 /M 2 ) 1 / 2 x 2 sinh 2 v? 1 
(1 - q 2 /M 2 ) 2 



q — > oo , x fixed , 



(7.17) 
(7.18) 

(7.19) 
(7.20) 



where 



tp = 2asinh 



if 2 sjx 2 -x 



, x — » X2 , — (? 2 fixed , (7.21) 



For the leading isospin moments we then have 



2-7TX 



iV-1 



R 2 (l- M 2 /q 2 ) 2 J S(z,q 2 ) 



dz- 



sinh 2 t sinh 2 6L 
— -I 



where 



4 



4M 2 sinh 2 — 
2 



5(z,(? 2 ) = y/(l-z)(l-z + 2z/R) 
Numerically one extracts the behavior 

for iV > 2 , 



ln(-q 2 /M 2 ) 



(7.22) 

(7.23) 

(7.24) 
(7.25) 

(7.26) 



in perfect agreement with (j4.1()j) . The dominant asymptotic piece comes from the 
u k + " contribution in (|7.23|) . and the dominant large q 2 behavior originates from the 
singularity near threshold. 
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7.2 The current structure functions 



We now turn to the current structure functions for which the non-trivial parts are 
more complicated to compute in the 1/n approximation than those for the spins. 
We have 



00 1 

w ab;cdef^^ ^ _ ^ _±_ _yab;cdef y^V\ _ yba;cdef y^V\ _|_ j£Cib;cdef -^y[r] ^ (7.27) 

and so for the various invariant isospin channels 

00 ^ 

w l {q 2 ,x)=Y J -A ] tf' x )' ( 7 - 28 ) 

r=0 

with 

W W = 2wf ] +2wt ] + wl +1 \ (7.29) 

w [ { ] = Wl r] -Wl r] , (7.30) 

«,M = -Wl r] -Wl r] . (7.31) 

At leading order one has the simple diagram with the 1-particle pole in the 
s-channel; this contributes to the structure functions only terms oc 5(1 — x). 

We denote the amputated two-current two-spin correlation function by y^.cde/ 

with coefficients in the 1/n expansion Tj^j, whose imaginary parts are proportional 

M ' \l] 

to the Wg . There are 3 types of diagrams contributing to T[. which involve one 

A propagator. One is the box diagram, another involves a vertex correction, and 

the last involves a spin self energy diagram. They can be conveniently combined 

together 8 to yield (in Euclidean space) 

T£l(p E ,q E )=8 J ^D,(p E ,q E ,k)D„(p E ,q E ,k)d(k + q E )B(p E -k)- 1 , (7.32) 
with 

D^p, q, k) = {k + q/2)^d(k) + (p + q/2)^d(p + q) . (7.33) 
So contracting over p, v 

E T i;i>^)=^(p)+8 J ^D(p,q,k)d(k + q)B(p-k)- 1 

(7.34) 

D(p, q, k) = (p + q/2) 2 d(p + qf + \ [d(p + q) + d(k)} 

+ [-(P - k) 2 +p 2 +pq- M 2 ]d(p + q)d(k) - J [q 2 + 4M 2 ]d(k) 2 , (7.35) 



one can conveniently use lattice UV regularization at intermediate stages 
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where F is a function of p 2 alone. Since we are only interested in the structure 
function only one term appearing in the cutting rule is relevant. Still in Euclidean 
space, and omitting pole terms oc 5(1 — x), the relevant term is 

E T i [1 i> E < « E ) ~ " 16 ^ m4 / A ^^^2 B (P E + 9 E ; M, m) (/+ + /_) , (7.36) 

where A is some ultraviolet cutoff and 

m = 2Mcosh«/2, (7.37) 
J± = 2£(p E ,g E ,A ;± (p E ,< ? E )), (7.38) 
[M 2 (p 2 +gp)+m 2 (g 2 + gp)] 

T (g ^ )2 V[(P + q) 2 + {m + Af) 2 ] [(p + q) 2 + (m - M) 2 ] .(7.39) 



So going over to Minkowski space 
W m = i MH ( W-ZM ) r i . (C ++ C_)si„h^ 



with kq defined through 



o (k 2 + vr 2 ) ^[W 2 - (m + M) 2 ] [W 2 - (m - M) 2 ] 

(7.40) 



ko W — M 

cosh — = (7.41) 

2 2M v 7 



and 



Cfc = _ 4M2 + 4 ^ + f_ _i ri + 2(- Pg + m 2 -2M 2 M^)l 
2(VF 2 -M 2 ) 2 VF 2 -M 2L v ; v ;J 

+d(K ± )-i(-g 2 + 4M 2 )d(ET ± ) 2 , (7.42) 
^ 2 _ I m2 ( m2 + + ™ 2 (? 2 + QP)] 

K ± ~ IP ^2 



T M — - (m + M) 2 ] [VF 2 - (m - M ) 2 ] . (7.43) 

Numerically for small x we find consistency 9 with the general result ()3.14j) : 

W 1] ~^2- a i(-9 2 ), 2 = 1,2, (7.44) 
x In x 



9 we did not yet confirm this analytically 
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with 



ai(V) 



1 

2^ 



sinh# 

-q 2 = 4M 2 sinh 2 -, 
2 



(2 + cosh 9) 



(7.45) 
(7.46) 



since a\(z) is the leading order 1/n contribution to the Adler function A\(z). In 
leading order 1/n the current vacuum 2-point function amplitude I\ = i\ with 10 



= — 



(z + 4M 2 



-B(r), r 2 = z. 



(7.47) 



We have not yet computed the 1/n contribution to the isospin structure 

[2] 

function wo- This also involves computing T 3 L J „ which is more complicated because 
it requires the evaluation of 2-loop graphs involving also A propagators. 



8 Summary and Conclusions 

In this paper we calculated the DIS structure functions (and their moments) in 
the family of the 2-dimensional O(n) nonlinear sigma models using standard field 
theory techniques available in any asymptotically free field theory model. In the 
special case of the 0(3) model we compared the results to the non-perturbative 
(bootstrap) determination of the same structure functions. 

The very good agreement between the results (in the intermediate energy range, 
where both the perturbative field theory results and the non-perturbative bootstrap 
results are expected to be valid) strongly indicates - once again - that standard field 
theory and the bootstrap define the same model. On the other hand this agreement 
provides some further, indirect proof for all the assumptions that are used in the 
derivations in both methods. 

The study of the structure functions has lead us to two interesting findings. 
First, we found that the isospin structure functions exhibit exact Bjorken scaling: 
for —q 2 ^oo the structure functions go to some non-trivially ^-dependent limits. 
Second, the exact small x asymptotics of the structure functions are shown to be 
different from the soft-Pomeron like fractional power behavior: the asymptotics of 
(x times) the structure function is logarithmic. 

In the first case we have obtained concrete results in the 0(3) model only but 
we think that our findings are more generally valid: it is probable that exact Bjorken 
scaling is due to the presence of the (infinitely many) higher spin conserved charges 
characteristic to integrable models. Also in the second case we believe that the small 

"'For z = V we have = ± + ^S^f, [^h] = ai« 
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x asymptotics we found here is valid in a more general setting. Whether something 
similar applies to QCD remains to be seen. 
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A 0(n) notations and identities 

The O(n) generators Q ab act in the defining (vector) representation as 

'gabyc] 



■ f ab yd 
lt cd V ■> 



where the generator matrix is 



f ab = 5 ac 5 bd _ 5 ad 5 bc 



This corresponds to the usual relation 



Q a ,V b 



in the n = 3 case with Q a = \e abc Q bc . 



ie abc yc 



The generator matrix in the 2-index tensor representation is 



ab j.ab rdv , ecu j_ab 



cd:uv b cu 



'Til. I U <*, 



dv 7 



and similarly for higher representations. 

Projector matrices in the 2-index tensor representation are 





ab:a'b' 



1 

ab:a'b' 



-8 ab 5 a ' b ', 
n 

_ ^jaa' gbb' ^ab' ^ba' 



l( 5 ™> § bb' +5 ab' § ba'\_l 5 ab § a'b' 
2 V J n 



They satisfy 



I 



-,ab:a'b' ^a'b'm'"^" 



= hi PI 



ab;a"b" 



,ab\a'b' 



8 aa' 5 W ^ 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



(A.5) 
(A.6) 
(A.7) 



(A.8) 
(A.9) 
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and 



with 



pf' ab = n (A. 10) 



nr 1 nr < n ~ l ) ^ (n-l)(n + 2) 

In the n = 3 case 717 = 21 + 1. 

6-index invariant tensors, antisymmetric in the last two index pairs are 

x ab;cdef = § ab ^cegdf _ gcf gde^ ^ 
yab;cdef _ ^ac^begdf _ ^ad^begcf _ ^ac^bf^de _|_ ^ad^bf^ce ^ (A. 13) 

(and Y ba ' cde f). The irreducible combinations are 

R ab;cdef = }_ x ab;cdef ? 

n 

R ab;cdef _ ^_ ^yab^cdef _ yba;cdef^j ^ 



D aft;ctfe/ 
it 2 



2 1 / \ 

^_j£db;cdef _ _ I y«b;ccie/ _|_ yba;cc(e/ \ ^ 
77- 2 V / 



which satisfy 



with 



R ab;cdef R a'b';cdef = ^ ^ p a6;a'6' (A ^ ?) 



f = 2(n-l), fi = f 2 = 4(n-2). (A.18) 

Note that in the n = 3 case the antisymmetric tensor representation coincides with 
the vector representation and indeed in this case 

}_^cxy ^duv R ab;xyuv _ pab.cd ^ ^g-j 

Finally we note the following identities. The generators of the vector represen- 
tation satisfy 

c^: = X>^f ;cd (A.20) 

1=0 
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with 

V = 2(n - 1) , V 1 = V 2 = n - 2 , (A.21) 

and similarly for the antisymmetric tensor representation 

2 

j.ax A>x j.ax ±bx \ rp j^ab;cdef i a 

l cd;uv l ef;uv l cd;uv z ef;vu ~ 1 I ' (A.ZZj 

1=0 

where 

T = 4(n - 2) , Ti = n - 2 , T 2 = 4 - n . (A.23) 

B Particle states, rapidity integrals 

The r-particle "in" states are characterized by the O(n) labels a±,. . . ,a r and the 
decreasing set of rapidities 0\ , . . . , 9 r and their normalization is 



m {a 1 ,9 1 ; . . . ; a r , 9 r \a\, 0\) . . . ; a r , 9 r ) m — 

(4vr) r 5 aia 'i . . . 5 ar <5(9[ - 9i) . . . 5{9' r - 9 r ) , 

corresponding to the completeness relation in the r-particle sector 



(B.l) 



(B.2) 



/o r -i 
d9 r \o-i, 0i] ... ] o r , 9 r ) m in (ai, 9i] . . . ] a r , 9 r 
-oo 

As usual, we introduce the set of positive rapidity differences 
ui = 9i — 02 , U2 = #2 - #3 , • • • u r -i = 9 r ^i - 9 r (B.3) 
and the r-particle invariant mass M^ r \u) with the definition 

r 

(u)e ±A = MJ2 e±f?l =E r ±P r = P r ° ± P} , (B.4) 



i=i 



where M is the mass of the O(n) particles. 

The r-particle rapidity integral can now be written 



/oo /-8i rO r -i roo r 

dOi / d9 2 ■ ■ ■ / d9 r = dA Vu {r) , (B.5) 
-oo J —oo J —oo J —oo J 
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where 



oo poo roc 

du\ / du2 • • • / dn r _i 
o Jo Jo 



The inverse transformation is 



where 



j3j — Uj r , j — 1, . . . j T 1 , 

ujk = Uj + Uj + i + • • • + Uk-i , 1 < j < k < r . 

fir = 0, 



and 



v\ = — in 



r-l 



1=1 



We note that 



u_L + v_ = in /i r , 



where [x T = (u) /M is the dimensionless invariant mass. 

C S-matrix asymptotics 

The Zamolodchikov O(n) S-matrix is |1] 

Sab-U ) = vi(0)5 ab 5 cd + a 2 (0)6 ac 6 bd + a 3 (9)5 ad 6 bc . 



where 



am 



-2m9 s 2 (9) 



a 2 (9) = (n-2)9 
a 3 (9) = -2vri 



i-K -9 (n- 2)9 - 2iti 

82 (e) 



(n - 2)9 - 2m ' 

S2(9) 



(n - 2)6* - 2m ' 
and the "isospin 2" phase shift is given by 



roo ^ 

s 2(#) = — exp {2i I — sin(0u;) K n (uj) 
Jo w 
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with 



Specially for n = 3 



e -7TU, + 



K n(«)= x _ w • (C6) 



Using the asymptotic formula 

J™ ^ sin(^) k{u) £ | fe(0) + M + O (1) , (C.8) 
we get the large 9 asymptotics of the S-matrix, which can be written as 

W) = 5 ac 5 M + j^)6 tgg + O (1) . (C.9) 

D Residue asymptotics 

For any (r + 2)-particle form factor in the O(n) model Smirnov's residue axiom (Hj 
can be written as 



2i,, _ A 



^l,..a r (/? + »7T + . . . , r ) ^ 7 (5 afe ^... ar (01, . ... ft 

Sba\...a r ;b\...b r a 

where 

5 , 6oi...ar;bi-.6 r a(^l) • • • j $r|/?) = SbanahiP ~ Ql)Sc 1 a 2 ;c 2 b 2 {f3 — #2) " " 
■ ■ ■ S Cr _ 2ar _ 1 - jCr _ 1 b r _ 1 (P — 9r—l)S Cr _ 1 a r ;ab r {(3 ~ @r) ■ 

If (3 is large we can use (|C.9|) to get 

rL 1 ...a r (P + ™ + e,f3,9 1 ,...,9 r )^ 



(D.l) 



(D.2) 



Ci •^...6j...a r (^i' • • • ' ^) • 



(n - 2)e/3 ^ 

So far the operator index A did not play any role. For the case of tensor 
operators, where A is an O(n) (multi) index, the form factors are invariant tensors 
and (|D.3|) can equivalently be written 



4-7T 

+ + e, /?, 0i, . . . , r ) - (n _ 2)g/g t& ^,.. ar (^1, • • • , B r ) , (D.4) 
where t°^ B is the O(n) generator in the appropriate representation. 
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E Operator basis 

In the operator product expansions we use a basis spanned by A$ ab , B^ ,l ^ ab and 
g(J,l)ab ^ w j iere ^ggg ^agig elements are hermitian local operators at (0, 0) and / = 0, 1 
or 2 tensor operators (in the index pair ab) under O(n). Under the action of the 
parity operator V they transform as 



VA^ ab V = A^ ab , vB^ ah V = Bl J ' l)ab . 

Their Lorentz spin can be read off the relations 



M,A^ )ab 



0, 



(J,l)ab 



= -iJB. 



(E.l) 

(E.2) 
(E.3) 
(E.4) 



where J is a positive integer and Ai is the Lorentz boost operator. Finally under 
the action of the (anti-linear) CPT operator 6 



QA^ ab Q = (-l) l A^ )ab , 



QB (J : l)ab e 



D l B ( : mab . 



(E.5) 
(E.6) 
(E.7) 



The one-particle matrix elements of the above operators are parametrized as 

(a,9\AV cd \b,e)=P^ cd A«\ (E.8) 



( a ,e\Bi J » cd \b,e) 



(a,6\B ( J ,l)cd \b,9} 



pab;cd g(J,l) 



iM , 



(E.9) 
(E.10) 



Note that we have considered operators with non-vanishing one-particle matrix 
elements only. 



F Notations and Conventions 

We will use the notation 

W 2 = -iW 



(F.l) 
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for any vector (and higher tensor) index. The light-cone components are 



and similarly 



W± = ~(W TW 1 ) = ~(iW 2 TW 1 ) (F.2) 



d± = \{d Tdi) = \{id 2 Td x ) . (F.3) 



We treat the two-dimensional Euclidean coordinates exceptionally since here we use 

y± = Tyi~W2, (F.4) 

which gives 

y+y- = -y 2 , (f.5) 

as opposed to the Euclidean square of vectors in (|F.2|) , which is given by 

(W 2 ) E = Wl + Wf = -AW+ W- . (F.6) 
Two-dimensional Fourier transformation is indicated by tilde: 

f(Q) = J d 2 ye lQ yf(y). (F.7) 

For functions depending on y 2 only we also define 

f(Q) = (Q 2 Y +1 (Jp) J J d 2 y e** f(y 2 )• (F.8) 

G Asymptotic expansions 

Assume that S(y) has an asymptotic expansion 

S(y) = f \ a - 1 + 0(\°) (G.l) 



in terms of the effective coupling A defined in (|5.66jl . Then S(Q) can be asymptoti- 
cally expanded as 

S(Q) = (2^)(-l) J J!(l - «r)/ A CT + O (A CT+1 ) (G.2) 
in terms of the effective coupling A defined by 

i + ^-lnA = ln]^. (G.3) 
A n - 2 A m 
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In the special case 



we have 



5(j/)=/ i + O(l), (G.4) 
A 



S(Q) = (2-7r)(-l) J J! f (1 + -1^ A) + 0(A 2 ) (G.5) 

i.e. in this case we also know the coefficient of the sub-leading term. Finally if 

S(y) = /o In A + 0(A) (G.6) 

then 

S(Q) = -(2tt)(-1) j J! /o A + 0(A 2 ) . (G.7) 

An alternative way of presenting the above results is as follows. If the derivative 
of the function W(y 2 ) has the asymptotic expansion 

y2 ^ W{y2) = aX<T + ( G - 8 ) 
then in Fourier space we have 

W(Q 2 ) = -4to(-1) j JIX 7 + O (A CT+1 ) . (G.9) 
In the special case a = if 

y2 ^ W{y2) = « + ^ + O (A 2 ) (G.10) 

then 

W(Q 2 ) = -4vr(-l) 3 J! {a + pX + O (A 2 ) } . (G.ll) 
H Reduced spin and current form factor squares 

The space of homogeneous symmetric polynomials in r variables 0j , i = 1 . . . r is 

(r) 

spanned by products of a k , 1 < A; < r 



r (r) 



E (H.l) 

l<il<---<i fc <r 



For the reduced spin and current structure functions we have 
for r = 2: 
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= 4 (aj - 3a 2 ) + 8vr 2 , 
= a? - 2a 2 - 2vr 2 , 
= o\ - 6a 2 + 2vr 2 . 



(H.2) 
(H.3) 
(H.4) 



for r = 3: 



.(3) 

Jo 



4 (-6a?a 3 + 2a 2 af + 19aia 2 a 3 - 6af - 9af) 
+ 4vr 2 (4a 4 - 21a 2 a 2 - 19o-icr 3 + 34a 2 J ) 
+8vr 4 (9a 2 - 22a 2 ) + 64vr 6 , 

-2a?a 3 + a 2 a| + 2a 1 a 2 cr 3 - 2a 2 + 9a| 
-2vr 2 (erf - 3a 2 a 2 + 10aia 3 - 4af) 
-2vr 4 (3a 2 - a 2 ) - 8vr 6 , 

6a? a 3 - a 2 a| - 38aia 2 a 3 + 6a| + 99a| 
-2vr 2 (a 4 - 9a 2 a 2 + 8aia 3 + 16a|) 
-2vr 4 (3a 2 - 17a 2 ) - 8vr 6 . 



(H.5) 




(H.6) 



(H.7) 



for r = 4: 



jj 4) = \<oo\o\o\ - 48a|a 2 - 48a?a| + 152ai a 2 a| - 72a| 

-(48a 2 a 2 - 144a| - 152a? a 2 a 3 + 476aia 2 a 3 - 56a 2 af - 52a 2 a|)a 4 

-(72a 4 - 52a 2 a 2 - 352a| - 128ai a 3 )a| - 640a| 

+47r 2 [8a 2 a| - 24af - 42a?a|a 3 + 131ai a|a 3 + 68a 4 a| 

-249a 2 a 2 a 2 + 73afa 2 + 55aiaf - (38a 4 a 2 - 279a 2 a\ + 528a| 

+ 175a?a 3 - 669aia 2 a 3 + 133a|)a 4 - (157a 2 + 32a 2 )af] 

+47r 4 [36a 4 ai - 184a 2 a\ + 260a 2 - 88afa 3 + 447a? a 2 a 3 

+253a 2 a| - 696aia|a 3 - 248a 2 a| 

+ (185a 4 - 1156a 2 a 2 + 2120a| - 740aia 3 )a 4 + 96a|] 

+47r 6 [32a? - 232a 4 a 2 + 780a 2 a\ - 992af - 532a? a 3 + 1269ai a 2 a 3 

+295a| + (1211a 2 - 3220a 2 )a 4 ] 

+ 16vr 8 [60a 4 - 278a 2 a 2 + 405a| - 212aia 3 + 434a 4 ] 

+ 1287r 10 [17a 2 - 37a 2 ] + 1280tt 12 , (H.8) 
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cr l cr 2 cr 3 — ^cf "! — ^cfcf + 2ai<720"! + 9(73 

+ (-2(7? + 4(72 + 2crf (7 2 (7 3 + 2a 1 aj(T 3 + 18(7? o| - 62fT 2 Cr|)(74 

+ (9cr? - 627?7 2 + 80cr|)cr| + 128cr! 

— 27T 2 [(7i<72 — 2(72 ~~ 37i7 2 73 + 3717^73 — 47?7 2 + 337?727 2 — 397 2 7 2 

+7cricr| + (107?7 2 - 457?7 2 + 78(72 - 15cr?(7 3 - 36cricr 2 cr 3 - 118(7|)(J4 
+(-14(7? + 2647 2 )7 2 ] 

+ 27T 4 [-3(7f(7| + 12(7?Crf - 16(72* + (7^(7 3 + (7?(7 2 (7 3 + 587? of 

+ 9(71(7^(73 - 174(7 2 Cr| + (15(7? - 202(7? (7 2 + 456o| - 256(71 (7 3 )(7 4 + 1184(71] 

-2tt 6 [4(7? - 24(7^(72 + 57a? a\ - 26o"| + 37(7? cr 3 

-198(7i(72(7 3 + 26of + (-272(7? + 1100(7 2 )(74] 
+ 87T 8 [-3(7f - 2(7?7 2 + 13(7| + 3(7173 + 126(7 4 ] 

-256tt 10 (72 + 128^ 12 , (H.9) 

= <Ae\ol - 60"|(7| - 6(7?(7| + 38(7i(72(7| - 99(7^ 

+ (-6(7?(7| + 36(72 + 38(7^72(73 - 242717^3 + 2(7?Cr| + 670(7 2 (7|)(74 

+(-99(7? + 670(7?7 2 - 6567| - 1600(7K7 3 )(7l + 3584a| 

-2VT 2 [-(7?(7| + 6(7| + 9<7?<72<7 3 - 55(71 Cr|(7 3 - 16(7? erf 
+ 99(7?(7 2 7! + 103(7 2 (7| - 257(7i7f 

+ (-8(7?cr 2 + 21(7?cr| + 78(7^ + 1257?7 3 - 7627i(7 2 (73 + 1664oJ )(7 4 
+(248(7? - 824(7 2 )(7l] 

+2TT 4 [3afaj - 32(7? <rf + 64a| - 17(7^(7 3 + 1657?7 2 7 3 

-280(7? of - 30371(7^(73 + 632(72(7f 

+ (557? - 284(7? (7 2 + 160(7^ + 632(7K7 3 )74 - 960(71] 

-2vr 6 [-4(7? + 56(7? (7 2 - 2017? a\ + 238cr| - 377?7 3 
-367i7 2 7 3 + 472cr| + (-827? + 2607 2 )7 4 ] 
+8vr 8 [37? - 467? 7 2 + 937^ + 35ctk7 3 + 587 4 ] 
+64tt 10 [7? - 87 2 ] + 128tt 12 . 

(H.10) 
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I Structure function moments 



log 10 (-<z 2 /M 2 ) 


M (1) 


M (2) 

lvl 0;2 


M ( ' A) 


M (4) 

iw 0;2 


M (5) 

iw 0;2 


1V1 Q;2 


1 


1.1434 


0.3675 


0.01199 


0.0001163 


7.7E-7 


5.13E-9 


2 


0.5766 


0.5516 


0.1065 


0.006140 


0.00015 


2.31E-6 


3 


0.3531 


0.5360 


0.2384 


0.04244 


0.0036 


0.000177 


4 


0.2305 


0.4591 


0.3200 


0.1076 


0.020 


0.002358 


5 


0.1592 


0.3762 


0.3456 


0.1716 


0.053 


0.01093 


6 


0.1154 


0.3051 


0.3370 


0.2168 


0.092 


0.02811 


7 


0.08700 


0.2484 


0.3121 


0.2410 


0.129 


0.05156 


8 


0.06777 


0.2041 


0.2818 


0.2490 


0.158 


0.07691 


9 


0.05419 


0.1696 


0.2513 


0.2457 


0.177 


0.1006 


10 


0.04427 


0.1425 


0.2229 


0.2363 


0.188 


0.1196 



(r ) 

Table 1: Values of moment Mq ;2 
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olU V / / 


1;1 


1;1 


1;1 


1;1 


1;1 


1;1 


1 
1 


n ^71 7 


3899 
u.oozz 


D 0904^ 


n 00007^ 


1 7F 


^ OF 8 


o 
z 


fl 988"} 




U. 1UOO 


n 00973 


^ RF ^ 


3 3F fi 


q 
o 


U. 1 f OD 


U.0000 


n 91 07 
U.Z1U 1 


u.uzuzy 


8 /IF ^ 


8 3F K 


4 


0.1153 


0.5127 


0.2879 


0.05878 


0.00286 


-0.00050 


5 


0.0796 


0.4545 


0.3301 


0.1072 


0.0139 


-0.00069 


6 


0.0577 


0.3986 


0.3451 


0.1526 


0.0342 


0.00162 


7 


0.04350 


0.3501 


0.3427 


0.1885 


0.0601 


0.00848 


8 


0.03389 


0.3096 


0.3307 


0.2133 


0.0873 


0.02014 


9 


0.02710 


0.2760 


0.3141 


0.2282 


0.1123 


0.03531 


10 


0.02214 


0.2481 


0.2959 


0.2354 


0.1333 


0.05215 



Table 2: Values of moment M{ r J 



logio(-9 2 / 


/M 2 ) 


^3 A/f (2fe-l) 


M (2k) 
l^k=l IV1 0;2 


^3 A/f (2fe-i) 
l^k=i lvl i-i 


v^3 M (2k) 
l^k=l 1V1 1;1 


1 




1.155 


0.3676 


0.5921 


0.3823 


2 




0.683 


0.5577 


0.3937 


0.5466 


3 




0.595 


0.5786 


0.3872 


0.5767 


4 




0.571 


0.5690 


0.4061 


0.5710 


5 




0.558 


0.5587 


0.4236 


0.5610 


6 




0.544 


0.5500 


0.4370 


0.5528 


7 




0.528 


0.5410 


0.4463 


0.5471 


8 




0.508 


0.5300 


0.4519 


0.5430 


9 




0.482 


0.5161 


0.4535 


0.5395 


10 




0.455 


0.4984 


0.4513 


0.5357 



Table 3: Values of sums of moments M^l and in the even and odd channels 
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J One-loop 2d integrals 

We start with the 1-loop Euclidean integral with 2 internal scalar propagators with 
masses mi, 7772: 

f°° d 2 q 1 
The integral can be done analytically to obtain 



B(k;mi,m2) = In 



k 2 + m 2 , + \ k 2 + m 



2ir^(k 2 +m 2 _) (k 2 +m 2 + ) { ^k 2 +m 2 + - ^k 2 + m 2 _ 



(J.2) 



where 

m± = mi ± m<i . (J-3) 
For the equal mass case mi = 771,2 = M we have 

B(k)^B(k;M,M) = 1 ^Vg+gg+Vg 

2iry/k?(k 2 + 4M 2 ) VA; 2 + 4M 2 - VP 

= 6(0) = - J . , for A; 2 = 4M 2 sinh 2 ^ . (J.5) 
v y 4vrM 2 sinh6' 2 v ; 

Note B{k) is analytic in k 2 with a cut from —00 to -4M 2 . Also B{k) / for all k 2 
and 

B(k) ~ for fc 2 -► 00 , (J.6) 

It can be represented by the dispersion relation 

-AM 2 



B 



00 

00 ^ 



= — / d/c— =— , (J. 9) 

2vr 7 fc 2 + 4M 2 cosh 2 f V ; 

where we have substituted z = — 4M 2 cosh 2 | and noted z±ie corresponds to setting 
9 = iTT ± k with k > 0: 

-j- [6(»7T + «) ~ 6(i7T - «)] = A ^ ~\ - . (J.10) 

2vr7 7 v /J 4vrM 2 sinh re v ; 
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The inverse of B satisfies a once subtracted dispersion relation 

,2 r -4M* Biz + ie)' 1 -B(z-ie)- 1 



7,2 p-4M' 

BW - 1 = B(0) -. + _/ d2 : 



Noting 



1 

2vr7 



we have 

B(k)- 1 = 4vrM 2 



b(m + k) b(iir — k) 



l + 2k 2 / dn 



r(z — k 2 ) 



4irM 2 sinh n 

K 2 + 7T 2 



sinh 2 | 



o (k 2 + tt 2 ) (fc 2 + 4M 2 cosh 2 §) 



J.l General 1 loop integrals ("cutting rule") 

We consider an arbitrary 1-loop integral {Ya=i h = ty- 



m 



2i -i 



where 



The result is simply 



where 



h — ] fcj ; (J>n — 0) • 
j'=i 



7 ( k ) = E ^ ( 7 5 + 7 tf ) BiUj-^mj) , 

i<j 



ti- 



ll [(« + 'r? + 

r=l,r^i,j 



21 "I | 



and the momenta are given by 



li 



2« 



(mf-m 2 ) i , 

"ft + *j) " - — 12 — Ll H =F 72- V4 + 4m i^' e • 



Z 2 - 



2 i2 , 2 2 



"Note (g± + Z;) 2 = and q% = qf. 
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